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' at non-zero temperature. 
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1 Introduction 



In the last few years AdS / CFT has been used to analyze systems in which condensed- 
matter-like phenomena occur. Among them are superconductivity and superfluidity 
(for a review see [1]), non- Fermi liquids [2l|3], metamagnetism [11[5], the quantum Hall 
effect [6HT3]. and more. In this note we provide another example, namely that of a 
striped phase of fermionic matter at non-zero density. 

Some of the aforementioned phenomena have been investigated in the bottom-up 
approach, in which one does not start with a well-defined boundary theory, but rather 
with a toy model in the bulk which seems to capture the essence of the phenomena. 
Here we take the top-down point of view and study the D3-D7' model, a particular 
( 2 -fl)- dimensional system consisting of a D7-brane intersecting a stack of D3-branes 
with i^ND = 6. We use the probe approximation, in which the D7-brane is treated 
as a probe in the near-horizon background of the D3-branes. The low-energy and 
weak-coupling spectrum has charged states that are purely fermionic. It therefore 
seems like a good place to start looking for interesting phenomena. 

The D3-D7' system was originally proposed in [T3] as a holographic model for 
strongly interacting fermions in 2+1 dimensions. In [Hill] this model was used to 
explore the phases of the system at non- zero temperature and charge density in the 
background of a magnetic field. In [111112] it was shown that the system exhibits a 
gapped phase similar to a quantum Hall state at fixed ratios of the magnetic field to 
the charge density. 

Here we analyze the system in the ungapped phase and without a magnetic field. 
This system resembles a Fermi liquid in that the resistivity at low temperature scales 
as T^. However, the heat capacity at low temperature does not follow the Fermi 
liquid formula. 

We analyze the four- dimensional bulk quasi-normal modes, which correspond to 
the poles in the bi-fermion two-point functions of the three-dimensional boundary 
theory, including the current-current correlators. In particular, we will demonstrate 
that the system becomes unstable to the formation of a striped phase when the 
ratio of the charge density to the squared temperature is above a critical value. This 
instability can be effectively described by a f/ (1) gauge field interacting with an axion 
field in a constant background electric field. 

Striped phases are known to exist in condensed matter systems. Examples include 
charge density waves and spin density waves (for reviews see [T^fTB] ). but there 
are more complicated situations as well. What is common to them is that some 
observable of the system is modulated in one direction and hence the name striped 
phase. Holographic striped phases in 3 + 1 dimensions were analyzed in pn - l20] . 

Note added : As this paper was being completed [21] appeared which finds a similar 
instability in a different (2+l)-dimensional holographic model. 
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2 An instability in Maxwell-axion theory 



In [T7] Nakamura, Ooguri, and Park demonstrated that five-dimensional Maxwell- 
Chern-Simons theory becomes unstable in the presence of a background electric field 
due to a non-zero momentum mode that becomes tachyonic. As we will now show, a 
similar instability can occur in a four- dimensional Maxwell-axion theory, which has 
the Lagrangian density 

L = + ^{dj<^)' - lm2$2 + ^e'-'^'^mjFKL • (1) 

The equations of motion are given by 

{a + m^)<!>-^ae'''''^ FijFkl = (2) 

diF^^ -2ae^'^^^di{^FKL) = 0. (3) 

Expanding around a background with F = F^^^ and $ = 0, the linearized equations 
for the fluctuations fu and ip are given by 

djf'-2ae''^%{^Fl^l) = 0. (5) 

Consider a constant background electric field in the 2;-direction Fq^^ = E. The equa- 
tions for the fluctuations become 

{d^d'' + did' + m^)ip-4:aEf^2 = (6) 
dj^' + d,r = (7) 
djf'^ + dir -AaEe'^d^if = 0, (8) 

where fi,^ take values in (0,3) and i,j take values in (1,2). Multiplying the third 
equation by ejkd'' and using the Bianchi identity gives another equation involving 
only /i2 and if, 

{d^d^ + did')h2 - AaEdid'^ = . (9) 
For plane wave solutions of the form 

(/.(x, t) = <^e-*(") , /i2(x, t) = f,,e"'i-'-^^) , (10) 

we obtain the dispersion relation 

cu^ = k^ + -m" ± -Vm4 + 64a2E2^2 _ 

Therefore, there is a tachyonic mode in the range < /c < y/lQa'^E'^ — ra^. For a 
background magnetic field, on the other hand, there is no instability. We will find a 
similar instability in the four-dimensional bulk theory of the D3-D7' system. This will 
be interpreted as an instability of the boundary theory at non-zero density towards 
the formation of a striped phase. 
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3 Review of the D3-D7' model 



A holographic model for fermions in 2+1 dimensions can be constructed using a D7- 
brane and D3-branes that intersect in 2+1 dimensions. Since i^ND = 6 in this 
case, the massless spectrum of the D3-D7 open strings consists only of the desired 
fermions. The separation between the D7-brane and D3-branes in the transverse 
direction corresponds to the fermion mass. We decouple the fermions from the closed 
string modes by focusing on the near-horizon limit of a large stack of D3-branes and 
treating the D7-brane as a probe. According to the usual gauge/gravity duality, the 
dynamics of the probe D7-brane in the near-horizon D3-brane background captures 
the physics of the (2+l)-dimensional fermions interacting with a strongly-coupled 
(3+l)-dimensional gauge theory. 

3.1 Background 

We begin with the near-horizon background of the non-extremal D3-branes: 

L~^ds% = r^-h{r)dt^ + dx^ + dy^ + dz^)+r-^(^-^ + r^dnl^ (12) 
Fs = iL'^ {r'^dt Adx Ady Adz Adr + dn^) , (13) 

where h{r) = 1 — r^/r^ and = y^AirgJ^a' . For convenience, we work in dimen- 
sionless coordinates, e.g., r = Vphys/ L. This background is dual to A/" = 4 SYM theory 
at a temperature T = rT/{'nL). We parameterize the five-sphere as an S"^ x fibered 
over an interval: 

d^l = #2^cos^z/'(rfr]5^V + sinV(rffi^^V 

{dQ^^^f = dej + siYi^e id(Pi, (14) 

where 0<?/'<7r/2,0<6'j<7r, and < 0j < 27r. As il) varies, the sizes of the two 
S'^'s change. At = one of the S'^'s shrinks to zero size, and ai ip = 7r/2 the other 
S"^ shrinks. The 5^ x 5^ at ^/^ = 7r/4 is the "equator". 

3.2 Probe 

The D7-brane extends in t, x, r and wraps the two two-spheres. The D7 embedding 
is usually taken to be characterized by ip{r) and z{r). However, excitations around 
this embedding are tachyonic. To cure this we turn on fluxes through the two two- 
spheres labeled by the parameters /i and /2. With the correct choice of /i and /2, one 
gets a stable embedding. We also consider a non-zero charge density, by including the 
time component of the worldvolume gauge field ao(^); and a background magnetic 



4 



field b. The D7-brane action has a DBI term given by 

Sdbi = -Ty j d^x -([ei{g^u + 2Txa'F^y) 

= J dr (4 cos^ i; + ff) (4 sin^ V + /I) 



X y (1 + r^hz'^ + r^h^'^ - a'o^) (^1 + , (15) 

and a CS term given by 

= -A^/1/2 /" rfrrV(r) + 2Ar /" rfr c(r)6ao(r) , (16) 



where A/" = 4:7r'^L^TjV2,i and c(r) = ■(/'(r) - | sin (4?/'(r)) - V'oo + ^ sin(4z/'oo)S Note 
that c(r), and therefore ^/'(r), plays the role of an axion in this model. Indeed, we 
will find a modulated instability of the type described in the previous section. 
The asymptotic behavior of the fields is given by 

i>{r) ~ ^/'oo + ■mr^+ - c^r^- (17) 
z{r) zo + ^ (18) 

ao(r) ~ /i--, (19) 
r 

where the boundary value ipoo and the exponents A± are fixed by the fluxes /i, f2'- 
{fl + 4 cos^ V'oo) sin^ = {fl + 4 sin^ V'oo) cos^ ^00 (20) 



A 3 1 , , + 16 cos6 ^00-12 cosVoo 



The parameters m and correspond to the "mass" and "condensate" of the fun- 
damental fermions, respectively, and yU and d to the chemical potential and charge 
density, respectively^] Throughout the paper we will choose the fluxes /i, /2 such that 
A+ = — 1 (and A_ = —2). We will vary fi and fix /2 to satisfy this condition. 



^For the MN embedding that describes the gapped phase, one must also add a boundary term 
Sboundary = 2A/'c(rmi„ )6ao (r,„i„ ) . However, we wiU only discuss BH embeddings, for which the 
boundary term vanishes. 

^The physical charge density is given by dphysicai = Sn^L'^a'Tid. 
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3.3 Embeddings 

There are two kinds of embeddings in general: "black hole" (BH) embeddings, in 
which the D7-brane crosses the horizon, and "Minkowski" (MN) embeddings, in which 
the D7-brane ends smoothly outside the horizon. The latter exist only for /i = or 
/2 = and for a fixed charge density to magnetic field ratio (filling fraction). Here 
we will only consider BH embeddings with a vanishing mag netic field 6 = oE 

For the fluctuation analysis, it is convenient to use the coordinate m = r-r/r. The 
boundary is now at m = 0, and the horizon is at m = 1. The embedding equations of 
motion now read 



u^rxziu) ~ z{u) 



OqW =aoM = d—- 

tt g{u) 

u'd^{g{u)^'{u)) = ^d^G{u), (22) 

2g{u) 

where for convenience we have defined a'o and z' and the prime now denotes differen- 
tiation with respect to u. We have also defined d = The functions G(u) and g(u) 
are defined slightly differently than in 



G{u) = (/2 + 4cos^^)(/2 + 4sin^V^) 



h d^u^ + G-flf^h 



The presence of the horizon imposes the boundary condition 



^'(1) = ---^^^^. (24) 

The mass associated with a given solution is given by 

m = mr^^ = i4>iu) — ijjoo) ■ (25) 

The BH embedding corresponds to the metallic phase of the model. In particular, 
it has a longitudinal conductivity at zero magnetic field given by [H] 



N. 
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27r2 

At low temperature and non-zero charge density 



x/c/2 + 4G(l). (26) 



o-xx ~ ^ , (27) 



^The fluctuations of the MN embeddings (at zero temperature) were studied in [12]. 
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which is the correct behavior for a Fermi hquid. Note, however, that this behavior 
is universal for any probe brane, and in particular holds also in cases where the 
fundamental degrees of freedom include bosons [22] • The heat capacity for a class of 
similar theories was computed in [23] and was found to behave at finite density and 
at low temperature as C„ ~ which is different from a regular Fermi liquid. 



4 Quasi-normal modes 

In this section we will study the fluctuations of all the fields around the BH embedding 
described in the previous section. Due to rotational invariance, we can restrict to 
fluctuations propagating only in the x-direction; we thus consider fluctuations of the 
form /(u)e~*'^*+*'^^. We rescale the fluctuations 



bz = rroz, dat,x,y = — — , oe.^ = — = ^ , (28) 



and define u = uj/tt and k = k/rx- This will remove tt from all the fluctuation 
equations. It is worth recalling here the relationship to the physical momenta Uphys = 
ttTCj and kp^yg = irTk. We also define the function 

A{u) = 1 + /i(m)mV(m)' + h{u)u-^z\u)^ - 4{u)^ . (29) 



4.1 Fluctuation equations 



Expanding the action to second order in fluctuations, we obtain a set of coupled, 
linearized equations for the fluctuations. The equation of motion for 6ip is given by 



h 



dlG - ^{d^Gf^ + y 5. log G]^ 5^ 



-u^du 



A 



(l + hu 



4 i::r2 -12 



z 



L 
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(1 + hu-^z") co'+{l + hu-'z'^ - 4^) hp 



log Gz 5z' 



gz'i,' 



h 



-Cj' + hk' 



log G6ci^ + u du 



dhv? 
A 



bz — u^du 



btp 

ip z bz 



Au^ 



k—^ip'bcx — ikAa^ sm^{2ip)bdy . (30) 
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The equation of motion for Sz is given by 



-Y^z'ip'Sip' - J-d^\ogGz'5ip 



+^ [- (1 + ^2 + (1 + _ -/2) ^pj 



A 



+ kdz'dcx — u^du 



dh 



z' 5d', 



(31) 



The equation of motion for 5dt is given by 

" 4 

where 



4-/2 



H{u) = ^ (1 + /in-V2 + /inY^) 5d[-%\ogG5^lj + ^^'5ij'-4^z'5z'. (33) 



z''' + hu^ip''^) 5ex—ikv?6dydu [2c{u)] = . 

(32) 

dhu^ I I dh _i I 



Ah' • ' ' 1 

The equation of motion for bd^ is given by 

4 

dku(z'Sz - u^^'S^) - (1 + hu-'^z'^ + huW 6e, 

h'' ^ 

gu^{—5e!^ + k5d[) — iuu^ddydu [2c(u)] = 0. 
Finally, the equation of motion for 5dy is given by 



Au' 



u 



+ -du 

00 



(34) 



ikAsm\2^)d'o6ij + ide^u^du [2c{u)] - ^ (l + hu-^z'^ + hu^^'^) u^6dy 

' 4 

+^APSdy - u^du [gu^5d!y\ = . 

There is also a constraint coming from the gauge choice a„ = 0: 

k 

-LoH + -gv^ (-54 + k5d[) = . 
4.2 Method of solution 



(35) 



(36) 



In this subsection we will briefly describe how the above equations of motion are solved 
to find quasi-normal modes. Although it is true that in some cases the equations of 
motion decouple from one another, we do not separately discuss those special cases. 
The method that we will describe here works for any number of coupled or decoupled 
differential equations. 
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First of all, we are interested in normalizable solutions with infalling boundary 
conditions for the different fields. By inspection, one finds that all the fields, except 
for ScLf^ will have the same infalling behavior as m — )■ 1: 

Silj,6z,6ex,Say ~ (1 — (^37) 
6at ~ (1 - u)^-'^ . (38) 

We separate out this leading singular behavior. For example, 5tp = {1 — u)~''^ S^reg, 
where 6ipreg is regular at the horizon. 

We are interested in finding those solutions for which all the fluctuations have a 
normalizable behavior near the AdS boundary. Generically choosing boundary condi- 
tions at the horizon and shooting towards m — ?• would be hopeless. However, one can 
utilize the so-called determinant method [2^125] , which was succesfully applied also in 

fr2\ . Here we have five field equations for five fields {Sipreg, ^^reg, S(^x,reg, So.y,reg, ^O.t,reg) 

to solve, but they are subject to a constraint. Imposing the constraint for this system 
is equivalent to imposing it on the boundary conditions at the horizon and satisfying 
all the other equations of motion. Thus, there are two alternative routes one can 
choose for taking the constraint into account: 1) make use of the constraint to solve 
for 5a[ in terms of all the other fields or 2) only impose the constraint on the horizon 
boundary conditions for the temporal gauge field. Both of the routes are equivalent, 
but we found it numerically faster to follow the second path. Therefore, we will not 
lose any equations but make sure to take the constraint into account. According to 
the determinant method, we choose a set of linearly independent boundary conditions 
at the horizon, say, 

{JVreg, 5^reg, 5ex,reg, ^Cly^reg} = {(1,1,1,1), (1, 1, 1, -1), (1, 1, - 1 , 1 ) , ( 1, - 1 , 1 , 1 ) } . 

(39) 

As mentioned above, the boundary condition for 6dt^reg at the horizon is set by the 
constraint: By expanding the constraint as a power series around the horizon, at 
zeroth order, one finds a condition for 6at^reg{^) in terms of others (this condition is 
equivalent to the relationship one would find from the 6at equation of motion ( 132|) ). 
The derivatives at the horizon are set by the equations of motion. For example, 
Sijj'j.^g{l) is found by expanding the 6ip equation of motion (l30|l around the horizon, 
and Sa[^^^g{l) is set by the boundary conditions and the derivatives of the other fields. 

Finally, for any given momentum we solve the set of differential equations four 
times, corresponding to the four different boundary conditions in fl39|) . The interesting 

^Notice that the different faU-off guarantees that the temporal component of the gauge field 
vanishes at the horizon, as should be the case when the thermal circle pinches off. 
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object to look at is the determinant 



det 



"^x,reg 
^^y,reg 



^reg 
'^'^reg 
"^x,reg 
^"'y,reg 



III 

reg 

'^'^reg 
"^x,reg 
y,reg 



'^^reg 
"^x,reg 



""'y,reg 



(40) 



at the AdS boundary. For a given k, one then begins to scan over the complex 
valued energy u until a zero of the determinant is found. Once this is the case, one 
concludes that a normalizable solution has been found. There is a linear combination 
of the boundary conditions giving the desired normalizable solution, for which all 
fluctuations vanish at the AdS boundary. 

In practice, we start with a limit of the parameters such that the equations decou- 
ple and consider separately the different fluctuations. The quasi-normal modes are 
identified as the values of (w, k) where the contribution to the determinant changes 
sign. The accuracy of these positions is therefore determined by the resolution of the 
scan, which in our case is at least 10~^. Away from this limit the determinant is 
complex in general, and we look for positions at which its absolute value is smaller 
than some "numerical zero". The "numerical zero" is chosen so as to maintain the 
same accuracy as in the decoupled limit. 



4.3 Hydrodynamical mode and zero sound 

Generically, the equations for the fluctuations are all coupled. However, there are 
certain limits of the parameters (d, m, fc, /i) where some of the equations decouple 
from the rest. This enables us to associate, even in the range of parameters where 
the equations are coupled, certain modes to particular fields at least by name. But 
the reader should not get confused; when the system is in the range of parameters 
where the equations are coupled, this naming may be meaningless. 

Let us start with the case where m = and d = 0. For this particular case 
ip' = and cXq = 0, so the equation for 6e^ decouples from the rest. For Se^ we find 
a hydrodynamical mode [i.e., Cj{k — )■ 0) — )■ 0) associated with charge conservation, 
which is purely imaginary and behaves as 

Cj = -iDP + . . . , A; ~ , (41) 

where D is the diffusion constant. For low k, this is the important mode since it 
has the lowest imaginary value, corresponding to the collective mode with the longest 
lifetime. However, as k is increased this mode meets another mode which initially 
started at u}{k = 0) ~ —2i. After the modes meet they become a pair of complex 
modes with both real and imaginary parts; the real parts have the same magnitude 
but opposite signs. We depict this in Fig. [1] (left panel). 
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Figure 1: Left: Quasi-normal modes for btx at c? = 0, m = and z/'oo = '?i"/4. 
Blue (solid) lines are the imaginary parts and red (dashed) lines are the real parts. 
Right: Diffusion constant clS 3i function of Notice the linear behavior for small 
temperatures, \. 

When c? 7^ 0, but still rn = 0, the longitudinal vector and the (5i scalar are 
coupled but the quasi-normal mode structure (if Cj is not too large) is quite similar to 
the zero-density case. There is a purely imaginary hydrodynamical mode that meets 
at some finite k another purely imaginary, but non-hydrodynamical, mode, and for 
larger k there are two complex modes. This transition from purely imaginary modes 
to complex modes can also be regarded as a transition from a hydrodynamical regime 
to a coUisionless regime (see [25] for a similar phenomenon in a (3+l)-dimensional 
boundary theory). The values of both \Cj\ and k at the meeting point decrease with 
increasing d\ see the left panel of Fig. O 
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Figure 2: Left: Quasi-normal mode for the coupled [Se^^ 8z) system at m = 0, = 10^ 
and i\}oo = 1^ Right: Imaginary part of the quasi-normal mode of the 5ex, 5z system, 
with m = 0, d = 5 and ipoo = 7r/4. 



The imaginary part of the complex mode for 7^ is shown on the right of 
Fig. El The imaginary part decreases as k increases and then increases for very large 
k towards zero. As d is increased, the minimum of Im Cj shifts to larger k. 

Indeed, this mode has the properties of the holographic zero sound mode [23l[26] . 
As T — 7- at fixed charge d (and therefore d — )■ 00), the purely diffusive mode disap- 
pears completely and the complex pair of modes start at A; — j- 0, with an approximate 
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dispersion relation 

Cj = vk- iak"^ + ... . (42) 

A calculation of for all /i, gives f ~ in the large d limit, in agreement with |[23] . 
In this model, in addition to the quantum attenuation at zero temperature, there 
is very large damping for low momentum at non-zero temperature, turning the zero 
sound mode to a purely dissipative mode. 

However, as we will see in the next subsection, these results can only be trusted 
up to some critical value of beyond which there is an instability towards a striped 
phase. 



4.4 Stripes 

In this section we show that at large enough d, an instability occurs. This instability 
occurs only at some non-zero kphys, similar to what was found in [T7] but now in a 
(2-|-l)-dimensional theory. This indicates that the Fermi liquid is unstable and that 
the true ground state of the (2-|-l)-dimensional theory will be a striped phase, similar 
to a charge density and spin wave. 




Figure 3: Lowest quasi- normal mode for 5dy (left) and dip (right), with d = 0, m = 
and ipoo = 7r/4. 



We start with c? = and m = 0. In this case the transverse gauge field fluctuation 
5dy and the embedding scalar fluctuation 5ip decouple. The lowest quasi-normal 
modes are portrayed in Fig. [3l At A; = the lowest quasi-normal mode has the same 
value of Cj as the second quasi-normal mode of the longitudinal vector (equations (IMl) 
and (l35|l are the same in this limit). 

For (i > and m = the modes 5dy and dip couple to one another. Above a critical 
value of d, the lowest negative imaginary quasi-normal mode crosses onto the upper 
half-plane, at some value of k, and becomes a positive imaginary frequency mode, 
i.e., an instability (Fig. HI left panel). The critical value is dc — 5.5. There is a range 
in k for which the frequency is positive imaginary, and at some larger k it crosses 
again back onto the lower part of the complex frequency plane. The dependence of 
this range on d is shown in Fig. H] (right panel). In this range the system is driven 
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towards a striped phase, where both the transverse current and the fermion bi-hnear 
are spatially modulated. The latter, in particular, describes a spin density wave. This 
can be seen, for example, from the fact that at finite density a spatially modulated 
field couples linearly to the magnetic field through the bulk CS term in (1161) . 



t 




Figure 4: Properties of the unstable mode with m = and tpoo = vr/4. Left: Lowest 
purely imaginary mode for the coupled {Say, 6ip) system. The lowest curve is for 
d = 5, and the upper curve is for (i = 6. Right: Range of instability; kmm (blue) and 
kmax (red) are plotted against d. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 5: Domains of instability for non-zero mass m and ipoo = vr/4 for d = 6 (smaller 
domain) and d = 10 (larger domain). 

The situation is similar when rh ^ 0. In this case all the fluctuations are coupled, 
and we still find an instability. As \rh\ increases, the critical value of d grows; see 
Fig. [5l Since now 6dt is also involved, one also obtains a charge density wave. As 
T — !■ the instability appears for any small, non-zero charge density. The physical 
momentum at which this happens is given by kp^ys — Things are qualitatively 

similar when ip^o 7^ vr/4. For instance when one decreases ipoo from the value 7r/4 
then dc decreases slightly. 
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5 Conclusion 



The ungapped (BH embedding) phase of the D3-D7' system has an instabihty to- 
wards a striped phase above a critical value of the ratio of the charge density to the 
temperature squared. For a non-zero mass the instability involves all the fields, and 
therefore one expects that the new phase will be such that all the physical quantities 
(normalizable modes) associated with these fields are spatially modulated. These in- 
clude charge, spin, and the transverse currents. In the dual bulk theory the instability 
is the result of an axion-like term in the action. On the other hand, the existence of 
this term is intimately related to the fermionic nature of the boundary theory. We 
believe that striped phases are a general property of such theories. Since the model 
should have a striped phase, one should be able to construct it as a solution of the 
equation of motion and study its properties. 
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